Direct Determination of the Structure of Single Biopolymer Molecules Using Nanopore Sequencing.

This review highlights operational principles, features, and modern aspects of the development of third-generation sequencing technology of biopolymers focusing on the nucleic acids analysis, namely the nanopore sequencing system. Basics of the method and technical solutions used for its realization are considered, from the first works showing the possibility of creation of these systems to the easy-to-handle procedure developed by Oxford Nanopore Technologies company. Moreover, this review focuses on applications, which were developed and realized using equipment developed by the Oxford Nanopore Technologies, including assembly of whole genomes, methagenomics, direct analysis of the presence of modified bases.

Rapid identification of SARS-CoV-2 strains via isothermal enzymatic recombinase amplification and nanopore sequencing.

Surveillance of the SARS-CoV-2 genome has become a crucial technique in the management of COVID-19, aiding the pandemic response and supporting effective public health interventions. Typically, whole-genomic sequencing is used along with PCR-based target enrichment techniques to identify SARS-CoV-2 variants, which is a complicated and time-consuming process that requires central laboratory facilities. Thus, there is an urgent need to develop rapid and cost-effective tools for precise on-site detection and identification of SARS-CoV-2 strains. In this study, we demonstrate the rapid diagnosis of COVID-19 and identification of SARS-CoV-2 variants by amplification and sequencing of the entire SARS-CoV-2 S gene using isothermal enzymatic recombinase amplification combined with the advanced Oxford nanopore sequencing technique. The entire procedure, from sampling to sequencing, takes less than 8 hours and can be performed with limited resources. The newly developed method has noteworthy implications for examining the transmission dynamics of the virus, detecting novel genetic variants, and assessing the effect of mutations on diagnostic approaches, antiviral treatments, and vaccines.

Illumina and Nanopore sequencing in culture-negative samples from suspected lower respiratory tract infection patients.

Objective: To evaluate the diagnostic value of metagenomic sequencing technology based on Illumina and Nanopore sequencing platforms for patients with suspected lower respiratory tract infection whose pathogen could not be identified by conventional microbiological tests. Methods: Patients admitted to the Respiratory and Critical Care Medicine in Shanghai Ruijin Hospital were retrospectively studied from August 2021 to March 2022. Alveolar lavage or sputum was retained in patients with clinically suspected lower respiratory tract infection who were negative in conventional tests. Bronchoalveolar lavage fluid (BALF) samples were obtained using bronchoscopy. Sputum samples were collected, while BALF samples were not available due to bronchoscopy contraindications. Samples collected from enrolled patients were simultaneously sent for metagenomic sequencing on both platforms. Results: Thirty-eight patients with suspected LRTI were enrolled in this study, consisting of 36 parts of alveolar lavage and 2 parts of sputum. According to the infection diagnosis, 31 patients were confirmed to be infected with pathogens, while 7 patients were diagnosed with non-infectious disease. With regard to the diagnosis of infectious diseases, the sensitivity and specificity of Illumina and Nanopore to diagnose infection in patients were 80.6% vs. 93.5% and 42.9 vs. 28.6%, respectively. In patients diagnosed with bacterial, Mycobacterium, and fungal infections, the positive rates of Illumina and Nanopore sequencer were 71.4% vs. 78.6%, 36.4% vs. 90.9%, and 50% vs. 62.5%, respectively. In terms of pathogen diagnosis, the sensitivity and specificity of pathogens detected by Illumina and Nanopore were 55.6% vs. 77.8% and 42.9% vs. 28.6%, respectively. Among the patients treated with antibiotics in the last 2 weeks, 61.1% (11/18) and 77.8% (14/18) cases of pathogens were accurately detected by Illumina and Nanopore, respectively, among which 8 cases were detected jointly. The consistency between Illumina and diagnosis was 63.9% (23/36), while the consistency between Nanopore and diagnosis was 83.3% (30/36). Between Illumina and Nanopore sequencing methods, the consistency ratio was 55% (22/42) based on pathogen diagnosis. Conclusion: Both platforms play a certain value in infection diagnosis and pathogen diagnosis of CMT-negative suspected LRTI patients, providing a theoretical basis for clinical accurate diagnosis and symptomatic treatment. The Nanopore platform demonstrated potential advantages in the identification of Mycobacterium and could further provide another powerful approach for patients with suspected Mycobacterium infection.

Nanopore adaptive sampling of a metagenomic sample derived from a human monkeypox case.

In 2022, a series of human monkeypox cases in multiple countries led to the largest and most widespread outbreak outside the known endemic areas. Setup of proper genomic surveillance is of utmost importance to control such outbreaks. To this end, we performed Nanopore (PromethION P24) and Illumina (NextSeq. 2000) Whole Genome Sequencing (WGS) of a monkeypox sample. Adaptive sampling was applied for in silico depletion of the human host genome, allowing for the enrichment of low abundance viral DNA without a priori knowledge of sample composition. Nanopore sequencing allowed for high viral genome coverage, tracking of sample composition during sequencing, strain determination, and preliminary assessment of mutational pattern. In addition to that, only Nanopore data allowed us to resolve the entire monkeypox virus genome, with respect to two structural variants belonging to the genes OPG015 and OPG208. These SVs in important host range genes seem stable throughout the outbreak and are frequently misassembled and/or misannotated due to the prevalence of short read sequencing or short read first assembly. Ideally, standalone standard Illumina sequencing should not be used for Monkeypox WGS and de novo assembly, since it will obfuscate the structure of the genome, which has an impact on the quality and completeness of the genomes deposited in public databases and thus possibly on the ability to evaluate the complete genetic reason for the host range change of monkeypox in the current pandemic.

Long-read sequencing of extrachromosomal circular DNA and genome assembly of a Solanum lycopersicum breeding line revealed active LTR retrotransposons originating from S. Peruvianum L. introgressions.

Transposable elements (TEs) are a major force in the evolution of plant genomes. Differences in the transposition activities and landscapes of TEs can vary substantially, even in closely related species. Interspecific hybridization, a widely employed technique in tomato breeding, results in the creation of novel combinations of TEs from distinct species. The implications of this process for TE transposition activity have not been studied in modern cultivars. In this study, we used nanopore sequencing of extrachromosomal circular DNA (eccDNA) and identified two highly active Ty1/Copia LTR retrotransposon families of tomato (Solanum lycopersicum), called Salsa and Ketchup. Elements of these families produce thousands of eccDNAs under controlled conditions and epigenetic stress. EccDNA sequence analysis revealed that the major parts of eccDNA produced by Ketchup and Salsa exhibited low similarity to the S. lycopersicum genomic sequence. To trace the origin of these TEs, whole-genome nanopore sequencing and de novo genome assembly were performed. We found that these TEs occurred in a tomato breeding line via interspecific introgression from S. peruvianum. Our findings collectively show that interspecific introgressions can contribute to both genetic and phenotypic diversity not only by introducing novel genetic variants, but also by importing active transposable elements from other species.

Detection of ribonucleotides embedded in DNA by Nanopore sequencing.

Ribonucleotides represent the most common non-canonical nucleotides found in eukaryotic genomes. The sources of chromosome-embedded ribonucleotides and the mechanisms by which unrepaired rNMPs trigger genome instability and human pathologies are not fully understood. The available sequencing technologies only allow to indirectly deduce the genomic location of rNMPs. Oxford Nanopore Technologies (ONT) may overcome such limitation, revealing the sites of rNMPs incorporation in genomic DNA directly from raw sequencing signals. We synthesized two types of DNA molecules containing rNMPs at known or random positions and we developed data analysis pipelines for DNA-embedded ribonucleotides detection by ONT. We report that ONT can identify all four ribonucleotides incorporated in DNA by capturing rNMPs-specific alterations in nucleotide alignment features, current intensity, and dwell time. We propose that ONT may be successfully employed to directly map rNMPs in genomic DNA and we suggest a strategy to build an ad hoc basecaller to analyse native genomes.

Streamlining remote nanopore data access with slow5curl.

BACKGROUND: As adoption of nanopore sequencing technology continues to advance, the need to maintain large volumes of raw current signal data for reanalysis with updated algorithms is a growing challenge. Here we introduce slow5curl, a software package designed to streamline nanopore data sharing, accessibility, and reanalysis. RESULTS: Slow5curl allows a user to fetch a specified read or group of reads from a raw nanopore dataset stored on a remote server, such as a public data repository, without downloading the entire file. Slow5curl uses an index to quickly fetch specific reads from a large dataset in SLOW5/BLOW5 format and highly parallelized data access requests to maximize download speeds. Using all public nanopore data from the Human Pangenome Reference Consortium (>22 TB), we demonstrate how slow5curl can be used to quickly fetch and reanalyze raw signal reads corresponding to a set of target genes from each individual in large cohort dataset (n = 91), minimizing the time, egress costs, and local storage requirements for their reanalysis. CONCLUSIONS: We provide slow5curl as a free, open-source package that will reduce frictions in data sharing for the nanopore community: https://github.com/BonsonW/slow5curl.

Advancing pathogen surveillance by nanopore sequencing and genotype characterization of Acheta domesticus densovirus in mass-reared house crickets.

Rapid and reliable detection of pathogens is crucial to complement the growing industry of mass-reared insects, in order to safeguard the insect colonies from outbreak of diseases, which may cause significant economic loss. Current diagnostic methods are mainly based on conventional PCR and microscopic examination, requiring prior knowledge of disease symptoms and are limited to identifying known pathogens. Here, we present a rapid nanopore-based metagenomics approach for detecting entomopathogens from the European house cricket (Acheta domesticus). In this study, the Acheta domesticus densovirus (AdDV) was detected from diseased individuals using solely Nanopore sequencing. Virus reads and genome assemblies were obtained within twenty-four hours after sequencing. Subsequently, due to the length of the Nanopore reads, it was possible to reconstruct significantly large parts or even the entire AdDV genome to conduct studies for genotype identification. Variant analysis indicated the presence of three AdDV genotypes within the same house cricket population, with association to the vital status of the diseased crickets. This contrast provided compelling evidence for the existence of non-lethal AdDV genotypes. These findings demonstrated nanopore-based metagenomics sequencing as a powerful addition to the diagnostic tool kit for routine pathogen surveillance and diagnosis in the insect rearing industry.

Snapshot study of canine distemper virus in Bangladesh with on-site PCR detection and nanopore sequencing.

Canine distemper virus (CDV) is a highly contagious virus that affects domestic and wild animals, causing severe illness with high mortality rates. Rapid monitoring and sequencing can provide valuable information about circulating CDV strains, which may foster effective vaccination strategies and the successful integration of these into conservation programs. During two site visits in Bangladesh in 2023, we tested a mobile, deployable genomic surveillance setup to explore the genetic diversity and phylogenetic patterns of locally circulating CDV strains. We collected and analysed 355 oral swab samples from stray dogs in Rajshahi and Chattogram cities, Bangladesh. CDV-specific real-time RT-PCR was performed to screen the samples. Out of the 355 samples, 7.4% (10/135) from Rajshahi city and 0.9% (2/220) from Chattogram city tested positive for CDV. We applied a real-time RT-PCR assay and a pan-genotype CDV-specific amplicon-based Nanopore sequencing technology to obtain the near-completes. Five near-complete genome sequences were generated, with phylogenetic relation to the India-1/Asia-5 lineage previously identified in India. This is the first study to provide genomic data on CDV in Bangladesh and the first demonstration of a mobile laboratory setup as a powerful tool in rapid genomic surveillance and risk assessment for CDV in low resource regions.

Comprehensive identification of pathogenic microbes and antimicrobial resistance genes in food products using nanopore sequencing-based metagenomics.

Foodborne pathogens, particularly antimicrobial-resistant (AMR) bacteria, remain a significant threat to global health. Given the limitations of conventional culture-based approaches, which are limited in scope and time-consuming, metagenomic sequencing of food products emerges as a promising solution. This method provides a fast and comprehensive way to detect the presence of pathogenic microbes and antimicrobial resistance genes (ARGs). Notably, nanopore long-read sequencing provides more accurate bacterial taxonomic classification in comparison to short-read sequencing. Here, we revealed the impact of food types and attributes (origin, retail place, and food processing methods) on microbial communities and the AMR profile using nanopore metagenomic sequencing. We analyzed a total of 260 food products, including raw meat, sashimi, and ready-to-eat (RTE) vegetables. Clostridium botulinum, Acinetobacter baumannii, and Vibrio parahaemolyticus were identified as the top three foodborne pathogens in raw meat and sashimi. Importantly, even with low pathogen abundance, higher percentages of samples containing carbapenem and cephalosporin resistance genes were identified in chicken and RTE vegetables, respectively. In parallel, our results demonstrated that fresh, peeled, and minced foods exhibited higher levels of pathogenic bacteria. In conclusion, this comprehensive study offers invaluable data that can contribute to food safety assessments and serve as a basis for quality indicators.

CRISPR-Cas9-guided amplification-free genomic diagnosis for familial hypercholesterolemia using nanopore sequencing.

Familial hypercholesterolemia is an inherited disorder that remains underdiagnosed. Conventional genetic testing methods such as next-generation sequencing (NGS) or target PCR are based on the amplification process. Due to the efficiency limits of polymerase and ligase enzymes, these methods usually target short regions and do not detect large mutations straightforwardly. This study combined the long-read nanopore sequencing and CRISPR-Cas9 system to sequence the target DNA molecules without amplification. We originally designed and optimized the CRISPR-RNA panel to target the low-density lipoprotein receptor gene (LDLR) and proprotein convertase subtilisin/kexin type 9 gene (PCSK9) from human genomic DNA followed by nanopore sequencing. The average coverages for LDLR and PCSK9 were 106× and 420×, versus 1.2× for the background genome. Among them, continuous reads were 52x and 307x, respectively, and spanned the entire length of LDLR and PCSK9. We identified pathogenic mutations in both coding and splicing donor regions in LDLR. We also detected an 11,029 bp large deletion in another case. Furthermore, using continuous long reads generated from the benchmark experiment, we demonstrated how a false-positive 670 bp deletion caused by PCR amplification errors was easily eliminated.

Nanopore sequencing with T2T-CHM13 for accurate detection and preventing the transmission of structural rearrangements in highly repetitive heterochromatin regions in human embryos.

BACKGROUND: Structural rearrangements in highly repetitive heterochromatin regions can result in miscarriage or foetal malformations; however, detecting and preventing the transmission of these rearrangements has been challenging. Recently, the completion of sequencing of the complete human genome (T2T-CHM13) has made it possible to accurately characterise structural rearrangements in these regions. We developed a method based on T2T-CHM13 and nanopore sequencing to detect and block structural rearrangements in highly repetitive heterochromatin sequences. METHODS: T2T-CHM13-based "Mapping Allele with Resolved Carrier Status" was performed for couples who carry structural rearrangements in heterochromatin regions. Using nanopore sequencing and the T2T-CHM13 reference genome, the precise breakpoints of inversions and translocations close to the centromere were detected and haplotypes were constructed using flanking single-nucleotide polymorphisms (SNPs). Haplotype linkage analysis was then performed by comparing consistent parental SNPs with embryonic SNPs to determine whether the embryos carried hereditary inversions or balanced translocations. Based on copy number variation and haplotype linkage analysis, we transplanted normal embryos, which were further verified by an amniotic fluid test. RESULTS: To validate this approach, we used nanopore sequencing of families with inversions and reciprocal translocations close to the centromere. Using the T2T-CHM13 reference genome, we accurately detected inversions and translocations in centromeres, constructed haplotypes and prevented the transmission of structural rearrangements in the offspring. CONCLUSIONS: This study represents the first successful application of T2T-CHM13 in human reproduction and provides a feasible protocol for detecting and preventing the transmission of structural rearrangements of heterochromatin in embryos.

Nanopore sequencing in distinguishing between wild-type and vaccine strains of Varicella-Zoster virus.

BACKGROUND: The introduction of varicella vaccines into routine pediatric immunization programs has led to a considerable reduction in varicella incidence. However, there have been reports of varicella, herpes zoster, and meningitis caused by the vaccine strain of varicella-zoster virus (VZV), raising concerns. Establishing the relationship between the wild-type and vaccine strains in VZV infections among previously vaccinated individuals is crucial. Differences in the single nucleotide polymorphisms (SNPs) among vaccine strains can be utilized to identify the strain. In this study, we employed nanopore sequencing to identify VZV strains and analyzed clinical samples. METHODS: We retrospectively examined vesicle and cerebrospinal fluid samples from patients with VZV infections. One sample each of the wild-type and vaccine strains, previously identified using allelic discrimination real-time PCR and direct sequencing, served as controls. Ten samples with undetermined VZV strains were included. After DNA extraction, a long PCR targeting the VZV ORF62 region was executed. Nanopore sequencing identified SNPs, allowing discrimination between the vaccine and wild-type strains. RESULTS: Nanopore sequencing confirmed SNPs at previously reported sites (105,705, 106,262, 107,136, and 107,252), aiding in distinguishing between wild-type and vaccine strains. Among the ten unknown samples, nine were characterized as wild strains and one as a vaccine strain. Even in samples with low VZV DNA levels, nanopore sequencing was effective in strain identification. CONCLUSION: This study validates that nanopore sequencing is a reliable method for differentiating between the wild-type and vaccine strains of VZV. Its ability to produce long-read sequences is remarkable, allowing simultaneous confirmation of known SNPs and the detection of new mutations. Nanopore sequencing can serve as a valuable tool for the swift and precise identification of wild-type and vaccine strains and has potential applications in future VZV surveillance.

Target enrichment improves culture-independent detection of Neisseria gonorrhoeae and antimicrobial resistance determinants direct from clinical samples with Nanopore sequencing.

Multi-drug-resistant Neisseria gonorrhoeae infection is a significant public health risk. Rapidly detecting N. gonorrhoeae and antimicrobial-resistant (AMR) determinants by metagenomic sequencing of urine is possible, although high levels of host DNA and overgrowth of contaminating species hamper sequencing and limit N. gonorrhoeae genome coverage. We performed Nanopore sequencing of nucleic acid amplification test-positive urine samples and culture-positive urethral swabs with and without probe-based target enrichment, using a custom SureSelect panel, to investigate whether selective enrichment of N. gonorrhoeae DNA improves detection of both species and AMR determinants. Probes were designed to cover the entire N. gonorrhoeae genome, with tenfold enrichment of probes covering selected AMR determinants. Multiplexing was tested in a subset of samples. The proportion of sequence bases classified as N. gonorrhoeae increased in all samples after enrichment, from a median (IQR) of 0.05 % (0.01-0.1 %) to 76 % (42-82 %), giving a corresponding median improvement in fold genome coverage of 365 times (112-720). Over 20-fold coverage, required for robust AMR determinant detection, was achieved in 13/15(87 %) samples, compared to 2/15(13 %) without enrichment. The four samples multiplexed together also achieved >20-fold genome coverage. Coverage of AMR determinants was sufficient to predict resistance conferred by changes in chromosomal genes, where present, and genome coverage also enabled phylogenetic relationships to be reconstructed. Probe-based target enrichment can improve N. gonorrhoeae genome coverage when sequencing DNA extracts directly from urine or urethral swabs, allowing for detection of AMR determinants. Additionally, multiplexing prior to enrichment provided enough genome coverage for AMR detection and reduces the costs associated with this method.

Evaluation of Automated Magnetic Bead-Based DNA Extraction for Detection of Short Tandem Repeat Expansions With Nanopore Sequencing.

BACKGROUND: Long-read technologies such as nanopore sequencing provide new opportunities to detect short tandem repeat expansions. Therefore, a DNA extraction method is necessary that minimizes DNA fragmentation and hence allows the identification of large repeat expansions. In this study, an automated magnetic bead-based DNA extraction method and the required EDTA blood storage conditions as well as DNA and sequencing quality were evaluated for their suitability for repeat expansion detection with nanopore sequencing. METHODS: DNA was extracted from EDTA blood, and subsequently, its concentration, purity, and integrity were assessed. DNA was then subjected to nanopore sequencing, and quality metrics of the obtained sequencing data were evaluated. RESULTS: DNA extracted from fresh EDTA blood as well as from cooled or frozen EDTA blood revealed high DNA integrity whereas storage at room temperature over 7 days had detrimental effects. After nanopore sequencing, the read length N50 values of approximately 9 kb were obtained, and based on adaptive sampling of samples with a known repeat expansion, repeat expansions up to 10 kb could be detected. CONCLUSION: The automated magnetic bead-based DNA extraction was sufficient to detect short tandem repeat expansions, omitting the need for high-molecular-weight DNA extraction methods. Therefore, DNA should be extracted either from fresh blood or from blood stored in cooled or frozen conditions. Consequently, this study may help other laboratories to evaluate their DNA extraction method regarding the suitability for detecting repeat expansions with nanopore sequencing.

Icarust, a real-time simulator for Oxford Nanopore adaptive sampling.

MOTIVATION: Oxford Nanopore Technologies (ONT) sequencers enable real-time generation of sequence data, which allows for concurrent analysis during a run. Adaptive sampling leverages this real-time capability in extremis, rejecting or accepting reads for sequencing based on assessment of the sequence from the start of each read. This functionality is provided by ONT's software, MinKNOW (Oxford Nanopore Technologies). Designing and developing software to take advantage of adaptive sampling can be costly in terms of sequencing consumables, using precious samples and preparing sequencing libraries. MinKNOW addresses this in part by allowing the replay of previously sequenced runs for testing. However, as we show, the sequencing output only partially changes in response to adaptive sampling instructions. Here we present Icarust, a tool enabling more accurate approximations of sequencing runs. Icarust recreates all the required endpoints of MinKNOW to perform adaptive sampling and writes output compatible with current base-callers and analysis pipelines. Icarust serves nanopore signal simulating a MinION or PromethION flow cell experiment from any reference genome using either R9 or R10 pore models. We show that simulating sequencing runs with Icarust provides a realistic testing and development environment for software exploiting the real-time nature of Nanopore sequencing. AVAILABILITY AND IMPLEMENTATION: All code is open source and freely available here-https://github.com/LooseLab/Icarust. Icarust is implemented in Rust, with a docker container also available. The data underlying this article will be shared on reasonable request to the corresponding author.

Long-term monitoring of SARS-CoV-2 variants in wastewater using a coordinated workflow of droplet digital PCR and nanopore sequencing.

Quantitative polymerase chain reaction (PCR) and genome sequencing are important methods for wastewater surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The reverse transcription-droplet digital PCR (RT-ddPCR) is a highly sensitive method for quantifying SARS-CoV-2 RNA in wastewater samples to track the trends of viral activity levels but cannot identify new variants. It also takes time to develop new PCR-based assays targeting variants of interest. Whole genome sequencing (WGS) can be used to monitor known and new SARS-CoV-2 variants, but it is generally not quantitative. Several short-read sequencing techniques can be expensive and might experience delayed turnaround times when outsourced due to inadequate in-house resources. Recently, a portable nanopore sequencing system offers an affordable and real-time method for sequencing SARS-CoV-2 variants in wastewater. This technology has the potential to enable swift response to disease outbreaks without relying on clinical sequencing results. In addressing concerns related to rapid turnaround time and accurate variant analysis, both RT-ddPCR and nanopore sequencing methods were employed to monitor the emergence of SARS-CoV-2 variants in wastewater. This surveillance was conducted at 23 sewer maintenance hole sites and five wastewater treatment plants in Michigan from 2020 to 2022. In 2020, the wastewater samples were dominated by the parental variants (20A, 20C and 20 G), followed by 20I (Alpha, B.1.1.7) in early 2021 and the Delta variant of concern (VOC) in late 2021. For the year 2022, Omicron variants dominated. Nanopore sequencing has the potential to validate suspected variant cases that were initially undetermined by RT-ddPCR assays. The concordance rate between nanopore sequencing and RT-ddPCR assays in identifying SARS-CoV-2 variants to the clade-level was 76.9%. Notably, instances of disagreement between the two methods were most prominent in the identification of the parental and Omicron variants. We also showed that sequencing wastewater samples with SARS-CoV-2 N gene concentrations of >104 GC/100 ml as measured by RT-ddPCR improve genome recovery and coverage depth using MinION device. RT-ddPCR was better at detecting key spike protein mutations A67V, del69-70, K417N, L452R, N501Y, N679K, and R408S (p-value <0.05) as compared to nanopore sequencing. It is suggested that RT-ddPCR and nanopore sequencing should be coordinated in wastewater surveillance where RT-ddPCR can be used as a preliminary quantification method and nanopore sequencing as the confirmatory method for the detection of variants or identification of new variants. The RT-ddPCR and nanopore sequencing methods reported here can be adopted as a reliable in-house analysis of SARS-CoV-2 in wastewater for rapid community level surveillance and public health response.

Benchmarking long-read aligners and SV callers for structural variation detection in Oxford nanopore sequencing data.

Structural variants (SVs) are one of the significant types of DNA mutations and are typically defined as larger-than-50-bp genomic alterations that include insertions, deletions, duplications, inversions, and translocations. These modifications can profoundly impact the phenotypic characteristics and contribute to disorders like cancer, response to treatment, and infections. Four long-read aligners and five SV callers have been evaluated using three Oxford Nanopore NGS human genome datasets in terms of precision, recall, and F1-score statistical metrics, depth of coverage, and speed of analysis. The best SV caller regarding recall, precision, and F1-score when matched with different aligners at different coverage levels tend to vary depending on the dataset and the specific SV types being analyzed. However, based on our findings, Sniffles and CuteSV tend to perform well across different aligners and coverage levels, followed by SVIM, PBSV, and SVDSS in the last place. The CuteSV caller has the highest average F1-score (82.51%) and recall (78.50%), and Sniffles has the highest average precision value (94.33%). Minimap2 as an aligner and Sniffles as an SV caller act as a strong base for the pipeline of SV calling because of their high speed and reasonable accomplishment. PBSV has a lower average F1-score, precision, and recall and may generate more false positives and overlook some actual SVs. Our results are valuable in the comprehensive evaluation of popular SV callers and aligners as they provide insight into the performance of several long-read aligners and SV callers and serve as a reference for researchers in selecting the most suitable tools for SV detection.

pathMap: a path-based mapping tool for long noisy reads with high sensitivity.

With the rapid development of single-molecule sequencing (SMS) technologies, the output read length is continuously increasing. Mapping such reads onto a reference genome is one of the most fundamental tasks in sequence analysis. Mapping sensitivity is becoming a major concern since high sensitivity can detect more aligned regions on the reference and obtain more aligned bases, which are useful for downstream analysis. In this study, we present pathMap, a novel k-mer graph-based mapper that is specifically designed for mapping SMS reads with high sensitivity. By viewing the alignment chain as a path containing as many anchors as possible in the matched k-mer graph, pathMap treats chaining as a path selection problem in the directed graph. pathMap iteratively searches the longest path in the remaining nodes; more candidate chains with high quality can be effectively detected and aligned. Compared to other state-of-the-art mapping methods such as minimap2 and Winnowmap2, experiment results on simulated and real-life datasets demonstrate that pathMap obtains the number of mapped chains at least 11.50% more than its closest competitor and increases the mapping sensitivity by 17.28% and 13.84% of bases over the next-best mapper for Pacific Biosciences and Oxford Nanopore sequencing data, respectively. In addition, pathMap is more robust to sequence errors and more sensitive to species- and strain-specific identification of pathogens using MinION reads.